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Abstract 


NÁVAR, J. Spatial and temporal hydro-climatic variability in 
Durango, Mexico. Water Technology and Sciences. Vol. V, No. 1, 
January-February, 2014, pp. 103-123. 


This report presents information on the hydro-climate 
temporal and spatial variability for the State of Durango, 
Mexico. Annual rainfall data fitted the Gamma probabilistic 
density functions. The Gumbel-Fisher-Tipett density function 
fitted maximum daily precipitation. The log-Pearson Type 
III density function projected maximum annual discharge 
and daily maximum peakflows. Evapotranspiration was 
estimated by the Thornthwaite methodology. More than 80 
climatic stations with data lasting longer than 30 years were 
used to fit probability density parameters. Geographical 
prediction equations; e.g. mean annual rainfall = a + fl 
(Latitude) +82 (Longitude) + f3 (Altitude) estimated 
the spatial variability of precipitation, temperature and 
evapotranspiration. Reconstructed and instrumental long- 
term (1860-2004) time series (precipitation and streamflow) 
data was temporally analyzed using ARIMA and linear 
and nonlinear regression techniques. Results showed mean 
annual temperature and evapotranspiration are spatially 
variable, diminishing with an increment in latitude and 
altitude. Mean annual precipitation on the other side 
increases with latitude and altitude but reduces with 
longitude. This trend was reversed for the Pacific Ocean 
physiographic region. Hydro-climatic data showed a lack 
of linear monotonic patterns indicating the time series are 
stationary in the first momentum. Oscillations of several 
time-scales: seasonal, inter-annual, 3-7 years; and 9-12- 
year were consistently noted on the modeled reconstructed 
and instrumental time series, spectral density analysis and 
correlograms for annual precipitation and streamflow data. 
Large-scale synoptic climate events such as the seasonal 
monsoon, the inter-annual variability, the El Niño Southern 
and the Pacific Decadal Oscillations appeared to explain 
temporal oscillations. Spatial variability was observed in 
hydro-climatic maps constructed using GIS techniques. This 
spatial and temporal information can be used in the decision 
making of water management resources at the local and State 
spatial scales; and in studies of productivity, forest wildfires, 
and floods of the State of Durango, Mexico. 


Keywords: water resources, monsoonal effect, El Niño- 
Southern Oscillation, Pacific Decadal Oscillation, climatic 
variability, increased water use, raster and vectorial maps. 


Resumen 


NÁVAR, J. Variabilidad hidroclimática espacial y temporal en 
Durango, México. Tecnología y Ciencias del Agua. Vol. V, 
núm. 1, enero-febrero de 2014, pp. 103-123. 


Este reporte presenta información sobre la variabilidad temporal 
hidro-climática y espacial para el estado de Durango, México. 
Datos de precipitación anual se ajustaron a las funciones de 
densidad de probabilidad gamma. La función de densidad de 
Gumbel- Fisher-Tipett proyectó la precipitación diaria máxima. La 
función de densidad log-Pearson Tipo II simuló la descarga anual 
máxima y los caudales máximos diarios. La evapotranspiración 
se estimó mediante la metodología de Thornthwnaite. Más de 
80 estaciones climatológicas con datos de más de 30 años se 
utilizaron para ajustar los parámetros de las funciones de densidad 
probabilística. Ecuaciones geográficas, por ejemplo, precipitación 
media anual = a + 1 (latitud) + f2 (longitud) + p3 (altitud) 
estimaron la variabilidad espacial de la precipitación, temperatura 
y evapotranspiración. Datos de series de tiempo reconstruidos 
con técnicas dendrocronológicas a largo plazo (1860-2004) de 
precipitaciones y caudales se analizaron mediante técnicas de 
regresión lineales y no lineales, y por medio de modelos ARIMA. Los 
resultados mostraron que la temperatura y la evapotranspiración 
promedio anual son variables que disminuyen con un incremento 
en la latitud y altitud. La precipitación media anual, por otro 
lado, aumenta con la latitud y la altitud, pero se reduce con la 
longitud. Esta tendencia se invirtió en la región fisiográfica del 
océano Pacífico. Los datos hidroclimáticos mostraron una falta de 
patrones monótonos lineales, lo que indica que las series de tiempo 
son estacionarias en el primer momento. Oscilaciones de varias 
escalas de tiempo: estacionales, interanuales, 3-7 años, y de 9 a 
12 años se observaron consistentemente en las series de tiempo 
reconstruidas e instrumentales, el análisis de densidad espectral 
y los correlogramas para los datos de la precipitación anual y del 
caudal. Eventos climáticos sinópticos de gran escala, tales como 
los monzones estacionales, la variabilidad interanual, El Niño- 
Oscilación del Sur y La Oscilación Decadal del Pacífico parecían 
explicar las variaciones temporales. Se observó variabilidad 
hidroclimática espacial en mapas construidos utilizando técnicas 
SIG. Esta información espacial y temporal puede ser utilizada en la 
toma de decisiones en la gestión de los recursos del agua a las escalas 
espaciales locales y estatales, y en los estudios de productividad, 
los incendios forestales y las inundaciones del estado de Durango, 
México. 


Palabras clave: recursos hídricos, el efecto del monzón, El Niño- 
Oscilación del Sur, Oscilación Decadal del Pacífico, variabilidad 
climática, aumento de consumo de agua, mapas tipo ráster y 
vectoriales. 
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Introduction 


The diminishing of water resources by an 
increasing population, agriculture and industry; 
the impairment by pollution; and the erratic 
but consistent shortages by drought spells are 
issues society is currently and will be certainly 
facing during the next 50 years (Shiklomanov 
and Rodda, 2003; Postel and Wolf, 2001; 
Vorosmarty et al., 2000; Návar-Cháidez, 2011; 
Návar, 2012). Other current problems are due 
to: mismanagement, expansion of urban areas, 
economic development and poor decisions on 
water issues (CSIS, 2005). Therefore, in most of 
the arid and semi arid world or in places with a 
well defined wet-dry season there is currently 
a deficit between supply and demand that is 
critical for development. 

In general, Mexico's northern regions and in 
particular the State of Durango lies within the 
boreal hemisphere belt of arid and semi arid 
lands. Rivers of the state are currently running 
dry in several stretches due to the prolonged 
drought spell of the 2010's, in addition to 
river flow allocations to irrigated agriculture, 
industry, and public sectors. Groundwater 
is being heavily depleted resulting from 
increasing water demands and unsustainable 
development. At the present, in the region of 
Durango four aquifers are classified as over- 
drained (FAO, 2005). Moreover, drought spells 
of several temporal frequencies are further 
increasing water scarcity in the region (Návar, 
2012); e.g. the present inter-annual drought 
episode that is striking northern Mexico is one 
of the most acute ever recorded over the last 90 
years (CNA, 2011). 

In spite of the long-term availability of 
fresh renewable per-capita water resources in 
the State of Durango, Mexico, they are highly 
variable in space and time. The state has a range 
of discharge values between 13,000 to 15,000 M 
m? (11 800 M nr as surface flow and between 
1 000 and 3 000 M nY as aquifer discharge) per 
1.5 M inhabitants, which translates to close 
to 10 000 m? per inhabitant per year (Návar, 
2008). This figure is well above the Mexican 


present mean, which is close to 4 500-5 000 m* 
per inhabitant per year. However, the eastern 
portion of the State is arid and semi-arid with 
per-capita water resources well below the 
national average, with a figure that compares 
with that of the poorest Middle East countries. 
Drought spells of different magnitudes and 
durations reduce per-capita availability by 
more than one order of magnitude. The upper 
ridges of the Sierra Madre Occidental mountain 
range, on the other hand, have the largest per- 
capita water availability. 

In arid, semi-arid, sub-tropical, and tropical 
lands with a well defined dry season, Postel 
(2000) and CSIS (2005) stressed the need for the 
implementation of sustainable management 
practices of water resources to cope with spatial 
and temporal chronic water scarcity. Despite of 
these recommendations, there is a consistent 
lack of information publically available for 
the 
hydrologic projects. Since spatial and temporal 


development and implementation of 
information on water resources are required to 
better understand potential inherent sources 
of variability, geographic information system, 
GIS, coupled with statistical analysis of hydro- 
climatic variables are necessary for the display, 
management and forecasting of hydro-climatic 
data. This information is useful to take the 
appropriate measures during the presence of 
drought spells in places with consistent water 
shortages. 

In this paper, hydro-climatic information 
was statistically analyzed and compiled into 
a GIS system for the sustainable management 
of water resources of the State of Durango, 
Mexico with the aim to manage water resources 
for the well being of present and future 
generations. The aim of this research was: a) to 
geographically analyze hydro-climate variables 
by fitting equations between precipitation, 
temperature, and  evpotranspiration and 
coordinates (latitude and longitude) as well 
as altitude above sea level; b) to fit several 
density distribution functions, c) to statistically 
analyze whether hydro-climate (precipitation 
and river discharge) data shows consistent 
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temporal trends and oscillations that could be 
related to climate change and/or variability; 
and d) to discuss sources of inherent hydro- 
climatic pulses and oscillations. The intrinsic 
null hypothesis that 
variables do not present spatial or temporal 
variability in the state of Durango, Mexico. The 
spatial variability was tested by developing 
geographical equations to predict hydro- 
climate variables. The temporal variability was 


was hydro-climate 


evaluated by fitting autoregressive integrated 
moving average, ARIMA, models as well as 
linear regression equations to statistically 
analyze by monotonic trends in the first 
momentum. 


Materials and methods 
Location of the study area 
The State of Durango is located in Mexico's 
north-central portion covering an area of 


approximately 12.3 M ha (figure 1). It spans 22” 
35' NL and 104” 50” WL; 24” 44” NL and 22” 58” 


WL; 26” 83” NL and 104” 27” WL, and 23” 52 
NL and 107” 21” WL. It neighbors the States of 
Chihuahua and Coahuila to the north and east; 
Coahuila and Zacatecas to the east, Zacatecas 
and Nayarit to the south and Sinaloa and 
Nayarit to the west. Four main physiographic 
regions characterize the State; a) the western 
Plains of the Pacific Ocean, b) the Sierra Madre 
Occidental mountain range, c) the central 
valleys of Durango and Chihuahua, and d) the 
Chihuahuan Desert. 


Methods 


Instrumental climatic data for over 80 stations 
was available for the period of 1940-2004. Of 
these, 65 climatic stations are placed within the 
State of Durango and the rest within the State 
of Sinaloa, Mexico. In addition, daily discharge 
data was available for over 30 gauging stations 
placed most of them within the State and 
several of them just immediately outside the 
State of Durango. 


Pacific 
Ocean 


United States of America 


Golf of 
Mexico 


Figure 1. Location of the State of Durango in México. 
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Annual rainfall depth 


Amnual rainfall data fited the gamma density 
function with the aim to estimate annual 
precipitation depth with probabilities of 20, 
50 and 80%. The probabilistic density function 
that describes the gamma equation is given by 
the following model (Haan, 2002): 


px(x) = Mr T(n) X, An) 0 


Where x is the random variable of interest; 
e is the exponential mathematical function, A; 
y are the shape and scale parameters of the 
pdf; T (n) is the gamma function and it has the 
following properties: 


F(n) = (n -1)! paran=1,2,3,... 
T(n + 1) = nT(n) paran>0 
T(n) = foca paran>0 

0 


r(1)=T (2) =1:1(1/2) = Vx 


In order to fit this probabilistic density 
function it is necessary to estimate the shape 
and scale parameter values. They are quite 
frequently evaluated using the maximum 
likelihood solution technique of Greenwood 
and Durand (1975) coupled with a correction 
factor with the use of the following equations: 


y =Inx-Inx 
1 =(500876 +.1648852y -.0544274y?)/ y 


FC =3n/n 


Where In is the natural logarithm, x is 
the random variable of interest, e.g. monthly 
precipitation; A and y are the shape and scale 
parameters of the pdf; FC = correction factor. 


Rainfall depth-duration-frequency curves 


The Gumbel-Fisher-Tippet density function 
calculated rainfall depth for a 24-hour duration 
for the frequencies of 1, 2, 5, 10, 15, 25, 50 y 100 
years. The pdf of the Extreme Value Type l is 
given by the next model (Haan, 2002): 


Where u and a. are parameters of the pdf; xis 
the random variable of interest, e.2., maximum 
annual rainfall for a 24-hour duration; exp is 
the exponential mathematical function. The 
parameters of this pdf are estimated as follows: 


6 


a. =-—Sx; where Sx is the Standard deviation 
T 


of the random variable. 
u = m, 0.5772a; where m, is the mean of the 
random variable. 

The average annual number of days with 
measurable precipitation was also calculated 
for each of approximately 80 climatic stations. 


Potential evapotranspiration 


Annual potential evapotranspiration was 
preliminarily estimated using the Thornthwaite 
methodology. The empirical model developed 
by  Thornthwaite to estimate potential 
evaporanspiration is given by the following 
equation, as follows (Dunne and Leopold, 


1978): 


Erro a 


Where Et = potential evapotranspiration 
(cm); T, = mean air temperature ("C); I= annual 


13 15 
heat index = Ne , and a = 0.49 + 0.0179 
5 


¡1 
1 -0.0000771 1? +0.000000675 FP. 
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Base flow and direct flow 


Daily base flow and direct flow were separated 
using the RORA software (Rutledge, 1998) 
for 33 gauging stations. Daily discharge data 
recorded in all gauging stations are reported in 
the BANDAS program. 


Maximum hydrologic values 


Daily peak flows were also estimated for 
several probabilities of occurrence using the 
Log Pearson Type III density function. The 
Log Pearson Type II pdf has been widely 
recommended for these purposes (USGS, 
2002). The pdf for this model is given by the 
next equation, as follows: 


B 


A -(y-e) "-exp[-A- (y -e)] 


r(6) 


p(y)= 


The equations that solve for the parameters 
of the Log Pearson pdf by the moments 
methodology are the following; 


EIN 2 e=u _B 
Oy Cul YA 
Where: 
Sn(x, 
py =2 a ES 


” Flo] 


E 
n-1:n-2:0*? 


y 


This set of equations is solved simultaneously 
with the aid of computer programs. 


Reconstructed annual precipitation and 
discharge data 


Annual precipitation and discharge data 
were related to the annual radial increment 


of trees used typically in dendrochronological 
analysis. The dendrochronologies of Banderas, 
Las Bayas and El Gato reported by González- 
Elizondo et al. (2005) were used to relate the 
radial increment of trees to annual precipitation 
and discharge. Hence, regression equations 
reconstructed climate and discharge parameters 
beyond instrumental records using the past 
radial increment of trees (Stahle ef al., 1999). 
In general hydrologic data was reconstructed 
back to 1860. That is, the reconstructed and 
instrumental time series dates from 1860 to 
2004. 


Testing for monotonic trends 


In this report, for the variables annual 
precipitation and discharge, it is assumed 
that linear monotonic patterns coupled with 
oscillations are indications of climate change 
and variability, respectively. That is, should 
annual precipitation or discharge be controlled 
by climate change, the time series data must 
not be stationary in the first momentum. 
Hence, a simple linear regression equation (e.g., 
P = q + Bt) tests for the time series stationarity 
by setting the null hypothesis that B = 0. 
Second degree polynomial equations, e.g., 
quadratic patterns may be better indicators of 
climate variability. The argument against the 
null hypothesis could be that any linear or 
quadratic pattern is part of longer time scale 
climate variability. However, at this time for 
instrumental and reconstructed records there 
is no other statistical technique available for 
isolating the effects of climate change and 
climate variability on hydro-climatic variables. 
ARIMA models with linear and quadratic 
tendencies fitted the smoothed and raw data 
series, respectively. In order to look more 
closely to long term tendencies, first and second 
degree polynomial equations were fitted to the 
raw reconstructed and instrumental data (1860- 
2000). In addition, spectral density analysis and 
partial auto-correlogram functions improved 
our understanding on cycles and patterns of 
precipitation and discharge to be associated 
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to climate variability, climate change for the 
management of natural resources. 
Geographical information system 

Hydrologic maps were developed using 
annual and monthly rainfall, temperature, 
evapotranspiration, rainfall depth for 24-hour 
duration for several frequencies by predicting 
the spatial distribution of these values using the 
latitude and longitude coordinates and altitude 
as independent variables. A set of prediction 
equations are preliminarily reported in this 
paper as the main input in future updating 
maps to spatially display this information. 


Results 
Geographical models 


The spatial hydro-climatic variability was 
observed in GIS maps. Examples of the 
geographic models developed to predict 
hydro-climatic variables are reported in Table 1. 
Prediction equations were divided in two sets 
since the Sierra Madre Occidental mountain 
range divides the state in two major climatic 
provinces; the Pacific Ocean and the Central 
Valleys including the Chihuahuan Desert, each 
one with a unique set of climatic gradients. 

Geographic maps at the state scale showed 
large variability, which was partially reduced 
by climatically dividing the state into two 
physiographic regions, where the upper ridges 
of the Sierra Madre Occidental mountain range 
are both the continental and the hydro-climatic 
divide. For the interior physiographic region, 
Et and temperature diminish with latitude 
and altitude but increase with longitude. The 
thermal gradientis 1.1 *C per each 1 000 meters, 
diminishing from the Chihuahuan Desert to the 
upper ridges of the SMW mountain range. All 
precipitation variables increase with latitude 
and altitude but diminish with longitude. 
That is, the interior northwestern upland 
watersheds are wetter in the upper ridges of 
the SMW mountain ranges. 


For the Pacific Ocean physiographic region, 
Et and temperature diminish with longitude 
and altitude but increases with latitude. The 
average thermal gradient is 6 “C per each 1000 
my diminishing from the Pacific Ocean coasts to 
the upper ridges of the Sierra Madre Occidental 
mountain ranges by close to 12 *C. All pre- 
cipitation variables increase with longitude 
and altitude but diminish with latitude. That 
is, the wettest range is located in upland, upper 
ridges of the southeastern watersheds of the 
State of Durango. 

Mean annual precipitation is 645, 583, and 
575 mm when evaluated by the arithmetic, 
Thiessen polygons (figure 2) and the isohyets 
methodologies, respectively. Summer mon- 
soonal precipitation is typical in the region 
since more than 65% of the annual rainfall is 
restricted to the period of July to September. The 
number of days with measurable precipitation 
is on the average 64 (+ 4) but it increases to 94 (+ 
12) in the SMW mountain range and diminishes 
to 46 (+ 9) in the Chihuahuan Desert (figure 3). 
The average number of rainy days for the state 
of Durango, Mexico is depicted in figure 3. 

Raw annual precipitation, smoothed moving 
average data and fitted ARIMA models are 
depicted in figures 4, 5 and 6. Although annual 
precipitation is monotonically diminishing 
in the Plains of the Pacific ocean, the Sierra 
Madre Occidental mountain range and the 
Central Valleys of Durango, with the exception 
of the Chihuahuan Desert (figures 4, 5, and 6), 
the first degree polynomial equations do not 
have statistical significance. That is, the time 
series are stationary in their first momentum. 
Second degree polynomial equations better 
fitted annual precipitation for the climatic 
stations Fco. 1. Madero, Guanaceví, Rosario, 
and Potrerillos. The tendencies are to reduce 
(Fco. I. Madero, Rosario, and Potrerillos) and 
to increase (Guanaceví) showing the lack of 
consistent climate change effects on annual 
precipitation and appearing non-significant 
monotonic changes are most likely due to 
longer time scale climatic variability. That is, 
climate variability appears to be the rule when 
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Table 1. Examples of geographical prediction equations of climatic variables as a function of coordinates and altitude for the 


State of Durango, Mexico. 


Western slopes of the Sierra Madre Occidental mountain range 


Climate parameter Equation R? Sx Cv 
Et -3 748.8559-98.53 134*1t+67.39 232*1g-0.0670*al 0.35 94 mm 

Temp M -70.13988-1.4669*1t+1.18664*1g-0.0011*al 0.50 0.99 “C 4 

Pt< 20% 13 543+74.31406*1t-138.88607*lg+ 0.00728*al 0.11 204 mm 40 

Pt< 50% 20 221+117.70681*1t-209.2767*1g+0.00990*al 0.23 193 mm 27 

Pt< 80% 29 700+189.97134*1t-311.80979*1g+0.01809*al 0.30 225 mm Da 

PtM 22 021+128.4657*1t-227.66175*lg+0.04216*al 0.30 193 mm 23 

No. of rainy days 1 550.1540+8.11432*It-15.81674*1g+0.00461*a1 0.42 12 No. 21 

PtRP1y 480.57920+4.45236*It-5.11338*lg+0.0001563*al 0.03 11 mm 25 

PtRP 2 y 797.87780-2.7107*1t-6.03389*1g-0.00384*al 0.52 8 mm 10 

PtRP 5 y 1 087.84307-9.23365*1t-6.88060*1g-0.00748*al 0.59 14 mm 12 

Pt RP 10 y 1 275.2190-13.60427*1t-7.38623*1g-0.00990*al 0.57 19 mm 14 

PtRP 25 y 1 515.60199-19.08744*It-8.0681*1g-0.01296*a1 0.55 26 mm 15 

Pt RP 50 y 1 691.26833-23.1804*It-8.54348*1g-0.01519*a1 0.54 31 mm 16 

Pt RP 100 y 1 869.4952-27.2171*It-9.05694*lg-0.0174*al 0.53 37 mm 17 


Eastern slopes of the Sierra Madre Occidental mountain range, The Central Valleys and Chihuahuan Desert 


Et 4 099.3558+44.06216*lt-36.7461*lg-0.2496*al 0.77 80 mm 10 
Temp M 92.6339+0.5124*It-0.7324*1g-0.0060*al 0.73 a 12 
Temp Mn 71.23053+0.3547*lt-0.5796*lg-0.0059*al 0.61 2.40C 35 
Temp Mx 112.8249+1.3058*1t-1.0218*lg-0.0060*al 0.73 2D HE 5 
Pt < 20% -3 120.5348-42.8920*1t+40.6139*1g+0.0844*al 0.45 84 mm 32 
Pt < 50% -3 620.793-56.4930*1t+48.9796*lg+0.1281*al 0.54 92 mm 24 
Pt < 80% -4 548.905-71.0450*1t+62.0633*lg+0.1637*al 0.62 100 mm 18 
PtM -8 580.6623-176.0496*lt+123.8427*1g+0.2033*al 0.73 162 mm 2 
No. of Rainy days -844.8399-14.2974*1t+11.7323*1g+0.0092*a1 0.74 11 No. 18 
PtRP 1 y -244.7799-4.42932*It+3.52068*Ig+ 0.00348*al 0.57 5mm 19 
PtRP 2 y -76.4611-4.0371*1t+1.92483*1g+0.00732*al 0.48 6mm 14 
RUREIS:y -237.39352-5.77008*I1t+4.04752*1g+0.00527*al 0.29 11 mm ie 
Pt RP 10 y -243.2048-6.4875*1t+4.3781*1g+0.00507*a1 0.21 15 mm 22 
PtRP 25 y -444.2397-9.5201*lt+7.0755*lg+0.00649*a1 0.27 18 mm 22 
Pt RP 50 y -487.1182-10.6039*1t+7.8261*lg+0.00662*al 0.23 22 mm 24 
Pt RP 100 y -725.195-11.7571*1t+10.3221*lg+0.01184*al 0.30 24 mm 24 


Where: Et = evapotranspiration (mm); Temp= temperature; Pt = precipitation (mm); M = mean; Mn = minimum; Mx = maximum, RP = return 


period or frequency (years); lt = latitude; lg = longitude; al = altitude; y = year. The units for latitude and longitude are degrees in UTM 


coordinates and altitude is in meters above sea level. 


explaining the quadratic trends since they are 
contrasting for these climatic stations. 

The ARIMA models had 
single auto-regressive parameter for 


consistently 
one 
annual 


reconstructed and instrumental 


precipitation for most climatic stations with 


the exception of those placed in or close to the 
Chihuahuan Desert (Cuencamé, San Juan de 
Gudalupe, Simón Bolívar, Francisco 1. Madero, 
Guanaceví, Tepehuanes) that required two 
autoregressive parameters to explain over 50% 
of the total annual precipitation variability. 
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Figure 2. The Thiessen methodology for estimating average annual precipitation for the State of Durango, Mexico. 
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Figure 3. The average number of rainy days per year for the State of Durango, Mexico. 
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Figure 4. Patterns of recorded instrumental (1940-2004) and reconstructed (1860-1940) annual precipitation data for climatic 
stations placed in the physiographic province of the Sierra Madre Occidental mountain range and the Central Valleys of 
Durango and Chihuahua of the State of Durango, Mexico. The probability of fitting a first and second degree polynomial 
equation is also depicted in figure labels. 
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Figure 5. Patterns of recorded (1940-2000) and reconstructed (1860-1940) annual precipitation data for three climatic stations 
located in the physiographic province of the Chihuahuan Desert of the State of Durango, Mexico. The probability of fitting a 
first and second degree polynomial equation is also depicted in figure labels. 


Annual precipitation is hence more variable 
in the Chihuahuan Desert and the Central 
Valleys of Durango and Chihuahua and least 
variable in the Sierra Madre Occidental moun- 
tain range and the Pacific Ocean physiographic 
regions. 

ARIMA models with asingle auto-regressive 
parameter are sufficient to predict the smoothed 
annual time series. Forecasts project cycles in 
a different time domain but spectral density 
analysis indicate the time series show peaks in 
the frequency bands of 2-yr; 3-6-yr; and 9-11yr 
cycles, respectively (figure 7). These cycles 
are consistent across physiographic regions 
indicating the similitude of large synoptic 


events controlling regional climate variability 
across physiographic regions. 

The climatic stations Francisco I. Made- 
ro and Guanaceví (East Leeside) and 
Potrerillos, and Rosario (West ward side) 
present statistical significance on the fitted 
second degree 


polynomial equations. 


However, none of the climatic stations 
analyzed present statistical significance on 
the first degree polinomial equation. That 
is, oscillations larger than 70 years appear 
to explain the quadratic equations with 
statistical significance. In USA, Knight et al. 
(2006) linked this long-term variability to 
the Atlantic Multidecadal Oscillation. The 
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Figure 6. Patterns of recorded (1940-2000) and reconstructed (1860-1940) annual precipitation data for four climatic stations 
placed in the physiographic province of the western slopes of the Sierra Madre Occidental mountain range of the State of 
Durango, Mexico. The probability of fitting a first and second degree polynomial equation is also depicted in figure labels. 


potential climate change monotonic effect episodes can be noted during the last century; 
on annual precipitation is missing from this a) in the 1900's, b) in the 1920's, c) in the 1950's, 
analysis, since the first momentum of the and d) in the 1990's; with an intermediate dry 
time series changes in a quadratic fashion. spell of a shorter duration in the 1980's. 

The ARIMA models for the reconstructed Arithmetic air temperature has a mean (+ 
time series and the instrumental recorded confidence intervals) value of 17.3 *C (+ 0.9 *C). 
annual precipitations depict four major drought It increases to 19 *C in the Chihuahuan Desert 
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Figure 7. Partial autocorrelation functions for annual discharge (top figures) for the gauging stations of: (a) San Pedro; (b) 
Acaponeta; and (c) Salomé Acosta; as well as for annual precipitation (bottom figures) for the climatic statiaons of: 
(a) Tepehuanes; (b) El Salto; (c) Cuencamé), and (d) Tamazula for reconstructed and instrumental time series (1860-2004) of the 
State of Durango, Mexico. 


and the Pacific Ocean physiographic regions tendencies are absent showing mean annual 
but diminishes to 13.2 *C in the Sierra Madre temperature time series are stationary as well 
Occidental mountain range (figure 8). in this physiographic region. That is, the effect 
The thermal gradient approximates to 0.11 of climate change on mean annual temperature 
“C and to 0.60 *C per every 100 m in altitude is missing. For the Pacific Ocean province, 
above sea level in the Pacific Ocean and the mean annual temperature is decreasing with 
Chihuahuan Desert to the upper ridges of statistical significance (p = 0.0001) although 
the Sierra Madre Occidental mountain range, this observation requires additional evidence 
respectively. The difference in both gradients is before it can be further discussed as one climatic 
explained by the local effect of the Pacific Ocean station was relocated back in the 1990's. For the 
on temperature of the western lowlands of the Chihuahuan Desert, mean annual maximum 
Sierra Madre Occidental mountain range. temperatureisstatistically increasing (p=0.009), 
The temporal variations for mean annual mean annual temperature is also increasing 
temperature for three physiographic regions (p = 0.001) but the mean minimum is not 
of the state of Durango, Mexico are reported in significantly changing (p = 0.99). 
figure 9. 
The mean annual temperatures are not River discharge 
transiently changing for the SMW province 
(maximum with p = 0.43, mean with p = 0.11, The annual discharge has mean and con- 
minimum with p = 0.83). That is, the monotonic fidence interval values of 583 MmY and 220 
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Figure 8. Mean annual air temperature for the State of Durango, Mexico. 


Mnv, respectively according to river flow data 
for 25 gauging stations placed in the State. 
The long-term variation is close to 50% of the 
annual mean discharge stressing the effect of 
droughts and the presence of hurricanes or 
cyclones. Itis clear that rivers draining into the 
Pacific Ocean produce the largest and rivers 
draining into the interior Chihuahuan Desert 
produce the smallest annual discharges (table 
zh 

Gauging stations placed in rivers draining 
into the Pacific Ocean produce the largest and 
gauging stations placed on rivers discharging 
into the interior basins on the Chihuahuan 
Desert produce the smallest daily peak flows 
for most probabilities (table 3). 

Hydrograph separation analysis using daily 
discharge data indicate base flow and direct 
flow account on the average for 36 and 64%, 
respectively, for 17 gauging stations (table 4). 
Low base flow figures are the result of shallow 


perturbed soils and erratic and intense rainfall 
events. Durango is a State with the largest 
and smallest rains since more than 65% of the 
annual rainfall is composed of small-depth 
rains that do not produce significant runoff 
and the remaining 35% is composed of rains 
that produce both aquifer recharge and surface 
runoff. 

Discharge data is not monotonically 
diminishing over time for most gauging 
stations of the State of Durango (figures 10 
and 11). The time series are stationary in the 
first momentum missing the effect of climate 
change on stream flow-producing variables. In 
fact, the tendency of annual discharge is highly 
erratic due to the presence-absence of runoff- 
promoting rains such as cyclones, hurricanes, 
tropical depressions and high un-seasonal 
precipitation events in the area that masks 
any potential climate-change control over this 
variable. Several human-induced practices that 
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Figure 9. Patterns of annual air temperature for three physiographic regions of the State of Durango, Mexico. 
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Table 2. Annual discharge for 25 gauging stations with several probabilities estimated by employing the Log Pearson Type II 
density function. 


Gauging Annual discharge (Mm?) with a probability of occurrence (P(x < X)) Er 

station 10% | 20% | 30% | 40% | 50% | 60% | 70% | 80% | 9% | 95% | 99% 
Acaponeta 547 688 818 948 1 093 1262 1 474 17 SS 2905| 4544| 1359 
Acatitan + 77 115 162 224 31 440 665 1171 1886| 4641 394 
Baluarte 615 85 965 SE) 11837 LOS 1875 2SZS, 3141 4092| 6935 1760 
Guatenipa 435 596 756 927 1129 1380 1715 2 232 3233| 4460| 8394 1611 
Huites ENS 175 2194| 2610| 3090| 3674| 4439 0599) 7793| 10437| 18758| 4196 
Ixpalino 516 691 859 1 036 1238 1 484 1 803 2 280 3160| 4186| 7234| 4196 
Tecusiapa 285 309 381 458 546 655 800| 1023| 1455| 1990| 3763 761 
Sta. Cruz 285 441 611 805 1 047 1364 1 808 2532| 4020| 5950| 12549 1669 
Las Habitas 688 897 1095 1300 1534 1816| 2181 2727| 3742| 4946| 8673 1972 
La Bayona 23 34 46 59 76 99 132 191 330 551 1 658 123 
Las Tortugas 54 Ze 105 185 167 210 268 359 536 756 1 467 237 
Tamazula 229 296 360 425 499 588 701 867 1168 1510| 2498 636 
La Huerta 250 353 457 SAL 708 879 1109 1 469 DAA 3034| 5795 1 067 
.. 396 519 639 766 913 1094 1333 1699 2 405 3271 6 089 1179 
El Palmito 151 242 349 484 670 942 1379 DOS 4444| 8235| 28120| 1015 
El Pino 4 6 7 9 11 14 19 28 50 87 285 10 
El Saltito ls 24 38 56 82 121 183 302 602| 1082| 3305 25 
Fdez. II 543 646 743 844 963 1114 1324 1671 2427| 3504| 7916| 1191 
N. Mendoza 6 S) 12 16 2 28 Sí) 53) 85 128 285 34 
Peña Águila 5 9 14 21 30 45 69 118 249 477 1705 93 
a 92 142 195 256 331 428 562 RT 1210 1757| 3548 550 
San Felipe 3 7 12 20 32 53 90 172 425 932 4311 96 
San Pedro 1076 1374 1650 193 2245 2617| 3089 IV 5 023 6448| 10644| 2773 
V. Guerrero 1 3 a 8 13 20 33 61 140 286 1135 31 
Sardinas 62 103 151 208 283 386 539 805| 1405| 2270| 579% 530 


control river discharge such as the construction 
of dams « reservoirs; land use changes and 
overgrazing practices; increasing withdrawal 
of water resources from aquifers and rivers, and 
other sources of inherent climate variability are 
variables operating over time as well. 

Fitted ARIMA models explain a little 
over 50% of the smoothed (t = 3 yr) annual 
discharge data and these stochastic models 
include  consistently two  auto-regressive 
parameters for most gauging stations showing 
annual discharge is more variable than annual 
precipitation for the state of Durango, Mexico. 


The first auto-regressive (Lag 1) parameter 


is positive demonstrating the potential two- 
year effect of annual variability of high rains 
and drought spells on annual discharge in the 
area. However, the second (Lag 2) is negative 
showing that annual discharge restores to its 
previous conditions after the presence of inter- 
annual drought spells or inter-annual high 
discharge rates. 

Annual discharge cycles of 2-3-year; 4-6- 
year, and 9-10-year are also noted on the time 
(1860-2004), 
density analysis and ARIMA models. Similar 


series correlograms,  spectral 
cycles were described earlier for annual 


precipitation. 
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Table 3. Daily peak flows estimated by the Log Pearson Type III density function for 11 different probabilities for 17 gauging 


stations placed on rivers of the State of Durango, Mexico. 


Daily peak flows (m* s*) with a probability of occurrence (%) (X < x) 


Gauging 
station 

99.5 99 95 90 80 50 20 5 2 1 0.5 0.2 0.10 0.05 
A. Melgar 872 697 439 323 251 176 142 127 123 122 121 120 119 119 
El Palmito 3738| 3285| 2370| 1747| 1274 634 278 112 67 47 33 21 16 12 
Fernández II 462 388 268 205 163 113 85 69 64 61 58 56 55 54 
La Huerta 2206 1974| 149| 1158 895 515 273 138 94 72 56 41 33 26 
N. Mendoza 194 123 47 24 13 6 5 2 2 2 il 1 1 1 
P. del Águila 2195 1 472 572 257 118 25 5 1 0 0 0 0 0 0 
R. Salcido 436 396 294 214 147 56 15 3 1 il 0 0 0 0 
Río el Oro 1 005 988 920 821 691 394 152 40 17 9 5 2 1 1 
S. Acosta 1257 15 909 725 566 320 157 70 44 Sil 23 15 11 8 
San Felipe 793 663 411 256 154 48 11 2 1 0 0 0 0 0 
V. Guerrero 178 154 107 76 54 25 10 4 2 il 1 1 0 0 
Acaponeta 10987| 7985| 4021| 2435| 1586 779 439 286 238 214 196 178 167 159 
Acatitan 11298| 8891| 4786| 2690| 1475 391 81 14 5 3 1 1 0 0 
Las Habitas | 20247| 17358| 10959| 6626| 3681 833 110 9 2 1 0 0 0 0 
San Pedro 6088| 5069| 3370| 2455| 1846| 1102 685 450 368 323 289 253 231 212 
Santa Cruz 4842| 3936| 2410| 1597| 1073 490 214 94 61 45 35 25 19 15 


Table 4. Direct and base flow for 17 gauging stations placed on major rivers of the State of Durango, Mexico. 


Total So Base Direct Base Direct 
Watershed Gauging station | Years | discharge PE flow flow flow flow 
(Mm?) (%) (Mm) 

Sinaloa La Huerta 28 1 094.5 561.1 26 74 284.6 809.9 
Nazas Agustín Melgar 29 1 179.4 502.0 23 77 271.3 908.1 
Pre-Sanped El Saltito 34 236.9 218.7 222 78 52 184.8 
Nazas Fernández II 31 1 190.6 505.7 23 77 273.8 916.8 
Nazas El Palmito 19 1 361.3 580.5 21 Ya 285.9 1075.5 
Pre-Sanped N. Mendoza 36 31.9 18.7 66 33 21.1 10.5 
Pre-Sanped P. del Águila 43 94.9 103.9 55 45 SL 42.7 
Pre-Sanped R. Salcido 54 58.9 51.2 53 47 31.2 27.7 
Nazas S. Acosta 25 500.2 275.4 25 5 5 5752 
Pre-Sanped V. Guerrero 38 30.7 27.91 58 42 17.8 12.9 
Sinaloa Santa Cruz 50 1 651.3 682.5 84 16 1387.1 264.2 
Sinaloa Tamazula 36 633.9 223.2 25 75 158.5 475.4 
Pre-Sanped El Bejuco 36 150.0 54.5 24 76 36.0 114.0 
Pre-Sanped Los Habitas 18 2 030.4 774.2 29 71 588.8 1 441.6 
Pre-Sanped San Pedro 50 2 760.9 1 003.3 26 74 717.8 2 043.1 
Pre-Sanped Acaponeta 50 1 358.9 514.5 26 74 353.3 1 005.6 
Río Culiacán Guatenipa II Sil 1 568.0 665.3 24 76 376.3 1 191.7 
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Figure 10. Patterns of recorded (1940-2000) and reconstructed (1860-1940) annual discharge data for six gauging stations 
along the Mezquital and Nazas watersheds located within the physiographic regions of the Sierra Madre Occidental and the 
Chihuahuan Desert of the State of Durango, Mexico. 
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Figure 11. Patterns of recorded (1940-2000) and reconstructed (1860-1940) annual discharge data for three major rivers of the 


Pacific Ocean lowlands in the States of Durango and Sinaloa, Mexico. 


Discussion 


Mexican climate spans from the dry conditions 
of the northwestern Sonoran Desert where 
annual average rainfall is less than 100 mm 
to the wet tropical climates that characterize 
the forest regions of southern Mexico, where 
average annual rainfall can reach over 2000 
mm (Livermann, 1999). The State of Durango 
lies right in the center of this hydro-climatic 
This gradient 
originates in latitudinal belts of atmospheric 
which shift 
westerlies bring precipitation to northern 


gradient. strong rainfall 


circulation, seasonally. The 


Mexico in winter (Cavazos and Hastenrath, 


1990). The highs associated 
with stable conditions and the inter-tropical 


subtropical 


convergence zone and trade winds deliver 
summer rainfall to the central and southern 
regions of the country. Other key influences 
on precipitation are the fall hurricanes on both 
the Caribbean and Pacific coasts, the summer 
monsoons in the north, and the summer high 
pressure that disrupts the flow of moist air 
and creates a period of dry conditions known 
as midsummer drought. The mountains and 
varied topography of Mexico dominates other 
climatic influences. 

Most of the State of Durango, Mexico with 
the exception of the eastern-most tip lies within 
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the region called the monsoon region, which 
stretches along the Sierra Madre Occidental 
in north western Mexico and follows the 
continental divide through eastern Arizona, 
most of New Mexico and extends into southern 
Colorado in the United States (Comrie and 
Glenn, 1998). This area is characterized by 
the strong mid to late summer precipitation 
maximum with considerable less precipitation 
during the rest of the year. The narrow southern 
boundary of the monsoon region most likely 
blends into the inter-tropical convergence zone, 
ITCZ, related rainfall of central and southern 
Mexico. Storms causing precipitation during 
the winter come from the west and progress 
north to south while summer precipitation 
tends to progress south to north as a true 
monsoon system (Dettinger et al., 1998; Higgins 
et al., 1999). However, the area controlled by 
the monsoon is highly erratic and partially 
responsible for the seasonal and potentially 
the inter-annual oscillations in precipitation 
and discharge in watersheds of the State of 
Durango, Mexico. 

Several climatic trends and oscillations 
contributing to the monsoonal variability are: 
the El Niño- Southern Oscillation, ENSO, the 
Pacific Decadal Oscillation, PDO, the North 
Atlantic Oscillation, NAO, and the Little Ice 
Age. Reconstructed precipitation for Durango, 
Mexico using dendrochronological analysis 
(1380-2000) showed the tropical rainfall index 
(IRD correlates well with instrumental and 
reconstructed winter precipitation (r = 0.49 
and 0.55, respectively), reflecting the strong 
ENSO modulation of cool season climate over 
northern Mexico (Cleaveland et al., 2003). 
ENSO and PDO indices were statistically 
associated with long term hydro-climatic 
variability over The Sierra Madre Occidental 
mountain range (Návar, 2012). Due to the lack 
of long term information, the NAO and Little 
Ice Ages have not been tele-connected with 
hydro-climate variability in the region. Long 
term studies stress that climate is more variable 
over northern Mexico due to all the sources 
and interactions of variability described above. 


Hence, northern and southern Mexico often 
show opposite climate signals. Lake sediment 
analysis from central Mexico showed that Late 
Pleistocene-Holocene wet climate conditions 
prevailed about 6 000 yr BP, with large 
oscillations between 6 000 and 5000 yr BP and 
a period of marked aridity in Central Mexico 
about 1 000 yr BP (Metcalfe et al., 2000). 

The El Niño / Southern Oscillation, the PDO, 
the NAO, and the Little Ice Ages have been 
previously tele-connected with the monsoon 
variability (Cavazos and Hastenrath, 1990). 
However, to physically elucidate in more detail 
the potential controlling mechanisms of the 
monsoonal variability, the role of the Hadley 
circulation and the function of the inter-tropical 
convergence zone in the shape and strength 
of the subtropical circulation over the region, 
among others must be further studied (Comrie 
and Glenn, 1998). Cavazos and Hastenrath 
(1990), and Woodhouse (1997) discussed also 
how the subtropical anticyclones and associated 
upper-level ridges centered over the region 
and adjacent Atlantic and Pacific oceans exert a 
dominant influence year-round, modulated by 
the role of the mid-latitude westerly circulation 
in winter and the North American monsoon in 
summer. 

The monotonic climate change effect is 
missing because of the stationarity of the 
time series, although modeling the effect of 
climate change on river flow and precipitation 
by the IPCC (2007) forecasts northern Mexico 
is supposedly to reduce discharge and 
precipitation by 0-25% with global warming. 
Villers-Ruiz and Trejo-Vázquez (1997) also 
modeled the expansion of dry climates and 
consequently discharge on most Mexican 
ecosystems by climate change. Mulholland et 
al. (1997) and the IPCC (2007), using general 
circulation models, predicted northern Mexico, 
in particular the lower Rio Bravo/Rio Grande 
belt, would receive 10% less rainfall with 
increasing temperatures by global warming. 
Further information in the near future would 
shed light on the potential effect of climate 
change on hydro-climate variability in the 
studied region. 
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Other causes of discharge variability are the 
spatio-temporal abstraction of water resources 
directly from rivers and aquifers for irrigation 
of an expanding agriculture, population, 
and industry of the State of Durango. The 
mismanagement of natural resources, in 
specific plant and soil resources leading to 
soil compaction and reduction of plant cover 
appears to play also a role in diminishing 
discharge over time since several springs are 
seeing reduced baseflows in the region. 


Conclusions 


Hydro-climate is highly variable in the State of 
Durango, Mexico since most of the State is found 
within the region of monsoonal activity. This 
climate is characterized by erratic oscillations 
that have been associated with seasonal, 
inter-annual, ENSO, and PDO variability and 
potentially by NAO and Little Ice Ages as well. 
The role of the Hadley circulation, the inter- 
tropical convergence zone, the migration of the 
subtropical anticyclones and associated upper- 
level ridges, modulated by the role of the mid- 
latitude westerly circulation appear to exert a 
strong influence on hydro-climate variability 
in winter and the North American monsoon 
in summer. To preliminarily cope with this 
large variability, the sustainable management 
of water resources is an alternate practice that 
can meet economic growth, social stability 
and environmental compatibility and must be 
implemented to comply with the well being of 
present and future generations. Other projects 
aiming to increase water availability in the 
region are welcomed. 
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